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PREFACE 


This  final  technical  report  summarizes  the  work,  performed  under 

Contract  Number  F44620- 75-C-0088  at  Harvard  University  during  the 

period  June  1,  1975  - September  30,  1977. 

Work  in  the  area  of  laser  induced  breakdown  and  multiphoton 

absorption  in  ultraviolet  and  infrared  laser  window  materials  was  carried 

out  during  the  entire  period,  with  N.  Bloembergen  as  principal 

investigator.  A brief  summary  is  given  in  Section  A,  with  ample 

references  to  previously  issued  reports  and  publications.  A detailed 

account  of  the  results  obtained  during  the  final  semiannual  reporting 

period,  April  1 - September  30,  1977  is  contained  in  Section  B. 

Support  for  work  in  the  area  of  superconducting  materials  was 

discontinued  in  September  1976.  A final  report  on  this  work  was 

distributed  at  that  time.  In  Section  C the  principal  investigator. 

Professor  M.  Tinkham,  provides  a summary. 

Support  for  work  in  the  area  of  fracture  mechanics  was  discontinued 

in  September  1976.  A final  report  on  this  work  was  distributed  at  that 

time.  In  Section  C the  principal  investigator.  Professor  J.  R.  Hutchinson, 

provides  a brief  summary  of  this  area. 

The  investigators  contributing  to  this  final  report  are: 

Professor  N.  Bloembergen,  Principal  Investigator 
Professor  J.  R.  Hutchinson,  Principal  Investigator 
Professor  M.  Tinkham,  Principal  I nvestigator 
Dr.  J.  H.  Bechtel,  Former  Research  Fellow 
Dr.  H.  Lotem,  Research  Fellow 

Dr.  W.  Lee  Smith,  Former  Research  Fellow  and  Graduate  Student 
Mr.  P.  Liu,  Graduate  Student.  t ^ * * 
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The  views  and  conclusions  contained  in  this  document  are  those  of 
the  authors  and  should  not  be  interpreted  as  necessarily  representing  the 
official  policies,  either  expressed  or  implied,  of  the  Defense  Advanced 
Research  Projects  Agency  or  the  U.  S.  Government, 
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A,  SUMMARY  OF  WORK  ON  LASER  INDUCED  BREAKDOWN  AND 

MULTIPHOTON  ABSORPTION  IN  UV  WINDOW  MATERIALS. 

Professor  N.  Bloembergen,  Principal  Investigator. 

The  i trinsic  mechanisms  of  avalanche  ionization  and  multi- 
photon absorption  in  optical  materials  transparent  in  the  ultraviolet 
have  been  investigated.  This  work  is  relevant  to  the  development  of 
high  power  excimer  lasers  operating  at  ultraviolet  wavelength,  to 
elucidate  the  causes  of  failure  of  UV  windows  and  coatings. 

It  has  been  shown  that  the  damage  threshold  of  the  intrinsic 
window  material,  excluding  absorbing  impurities,  increases  by  a factor 
two  or  three  in  going  from  the  infrared  to  the  visible,  then  reaches  a 
maximum  and  decreases  in  the  ultraviolet  region.  This  qualitative 
trend  can  be  explained  by  the  creation  of  an  absorbing  carrier  plasma. 

The  threshold  at  lower  frequencies  is  mostly  determined  by  avalanche 
ionization,  at  higher  frequencies  by  multiphoton  absorption. 

The  breakdown  threshold  data  have  been  discussed  extensively 

12.  . . 3-5 

in  previous  technical  reports  * and  in  several  publications. 

This  material  will  not  be  reviewed  again,  as  Technical  Report 
No.  665  "Dielectric  Breakdown  Induced  by  Picosecond  Laser  Pulses" 
received  a wide  distribution. 

The  techniques  of  the  measurement  of  two-photon  absorption  cross 
sections  and  preliminary  results  have  also  been  discussed  in  the 
previous  semiannual  technical  reports,  as  well  as  in  several 
publications.  Extensive  results  on  two-photon  absorption  cross 

sections  in  a large  number  of  UV  window  materials,  as  well  as  in 
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several  harmonic  generating  crystals,  including  KDP,  ADP  and  others, 
have  been  obtained  more  recently.  These  are  presented  in  detail  in 
Section  B of  this  report.  This  same  material  has  also  been  submitted 

for  publication  in  the  Physical  Review. 

The  results  show  that  the  two-photon  absorption  cross  sections 
are  smaller  than  10'^  cm/W,  if  the  band  gap  of  the  material  is  larger 
than  twice  the  photon  energy.  It  may  be  concluded  that  intrinsic  multi- 
photon absorption  is  not  responsible  for  the  damage  occurring  at  present 
in  high  power  UV  laser  systems.  The  damage  in  coatings  and  windows 
in  these  systems  must  be  caused  by  one-photon  absorption  due  to 
impurities  and  color  centers. 

The  intrinsic  two-photon  absorption  cross  section  in  transparent 

materials  with  a bandgap  between  tio)  and  Z/ito,  lies  in  the  range  of  3x  10 

to  3 X 10"^  cm/MW.  It  is  fortunate  that  sapphire,  quartz  and  the 

harmonic  generating  crystals  have  cross-sections  at  the  lower  end  of 

this  range,  so  that  they  are  still  quite  useful  in  the  frequency  range 

Zhm>  E > fiti).  Further  details  are  given  in  the  next  section, 

gap 

The  following  publications  and  reports  were  supported  in  part  by 
this  contract. 

1.  W.  L.  Smith,  J.  H.  Bechtel  and  N.  Bloembergen,  Dielectric 
Breakdown  Induced  by  Picosecond  Laser  Pulses,  Technical  Report 
No.  665,  Harvard  University,  October,  1976. 

2.  W.  L.  Smith,  J.  H.  Bechb^l  and  N.  Bloembergen,  "Picosecond 
Breakdown  Studies:  Threshold  and  Nonlinear  Index  Measurements 
and  Damage  Morphology,  " in  National  Bureau  of  Standards,  Special 
Publication  No.  435,  edited  by  A.  J.  Glass  and  A.  Guenther, 

p.  321,  1975. 

3.  W.  L.  Smith,  J.  H.  Bechtel  and  N.  Bloembergen,  "Picosecond 
laser  induced  damage  morphology,  " Optics  Comm.  _18,  592,  1976. 


-5- 


W.  L.  Smith,  J.  H.  Bechtel  and  N.  Bloembergen,  "Picosecond 
laser-induced  breakdown  at  5321  and  3547  A:  Observation  of 
frequency  dependent  behavior,"  Phys.  Rev.  Bl5,  4039,  1977. 

W.  L.  Smith,  P.  Liu  and  N.  Bloembergen,  "Superbroadening  in 
HoO  and  D2O  by  self-focused  picosecond  pulses  from  a YAfG:Nd 
laser,"  Phys.  Rev.  A15,  2396,  1977. 

J.  H.  Bechtel  and  W.  L.  Smith,  "Pulse  width  fluctuations  in  mode- 
locked  lasers,  " Phys.  Lett.  55A,  203,  1975. 

J.  H.  Bechtel  and  W.  L.  Smith,  "Two-photon  absorption  in  semi- 
conductors with  picosecond  laser  pulses,  " Phys.  Rev.  3515, 

1976. 

W.  L.  Smith  and  J.  H.  Bechtel,  "A  simple  technique  for  individual 
picosecond  laser  pulse  duration  measurements,  " J.  App.  Phys. 
1065,  1976. 

J.  H.  Bechtel  and  W.  L.  Smith,  "Two-photon  fluorescence  in 
attenuating  media,"  J.  App.  Phys.  46,  5055,  1976. 

H.  Lotem  and  Cid  B.  de  Araujo,  "Absolute  determination  of  the 
two-photon  absorption  coefficient  relative  to  the  inverse  Raman 
cross  section,"  Phys.  Rev.  B16,  1711,  1977. 

E.  Yablonovitch  and  H.  S.  Kwok,  "Electrical  triggering  of  an 
optical  breakdown  plasma  with  subnanosecond  jitter,  Appl.  Phys. 
Lett.  583,  1975. 
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B.  ABSOLUTE  TWO- PHOTON  ABSORPTION  COEFFICIENTS  AT 
355  AND  266  nm. 

Professor  N.  Bloembergen,  Principal  Investigator. 


The  absolute  two-photon  absorption  coefficients  of  u.  v,  trans- 
mitting materials  have  been  measured  using  well-calibrated  single 
picosecond  pulses,  at  the  third  and  fourth  harmonic  of  a mode  locked 
Nd:YAG  laser  system.  Two-photon  absorption  coefficients  of  the  order 
of  lO”^  cm/MW  were  measured  for  alkali-halides,  and  10  ^ cm/MW 
for  harmonic  generating  crystals.  In  materials  with  band  gap  greater 
than  Zh(X),  no  nonlinear  absorption  could  be  observed.  Calculations 
based  on  the  Keldysh  theory  with  one  adjustable  parameter  agree 
quite  well  with  the  experimentally  observed  dispersions  in  the  two- 
photon  absorption  coefficients.  The  effect  of  multi-photon  absorption 
on  the  conversion  efficiency  of  harmonic  generating  crystals  and  on 
the  breakdown  of  u.  v.  window  materials  is  also  discussed. 


I HOBE 
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I.  Introduction 

The  first  theoretical  study  on  two- photon  absorption  (TPA)  process 
was  carried  out  by  Maria  Goppert-Mayer  in  1931.  With  the  advent  of 
the  ruby  lasers  Kaiser  and  Garrett^  observed  TPA  at  optical  frequencies 
in  1961.  Several  techniques  have  subsequently  been  used  to  find  TPA 
coefficients,  for  example,  by  measuring  the  intensity  dependent  trans- 
mission through  two- photon  excitation.  Several  comprehensive  reviews 

ui  3-5 

of  TPA  are  available. 

Absolute  TPA  coefficients  can  only  be  obtained  directly  by  well- 

calibrated  lasers.  However,  relative  measurements  can  be  done  quite 

simply  by  the  two- channel  technique.  ^ The  same  technique  may  also  be 

used  to  obtain  absolute  TPA  coefficients  by  utilizing  known  Raman  cross 

sections.  Since  the  TPA  coefficient  is  related  to  the  imaginary  part  of 

the  third  order  nonlinear  susceptibility,  parametric  mixing  experiments 

such  as  three -wave  mixing  which  are  capable  of  measuring  the  real  and 

imaginary  part  of  the  third-order  nonlinear  susceptibility  have  been  used 

8 

to  find  TPA  cross  section. 

TPA  studies  have  both  scientific  and  technological  interest. 
Because  the  TPA  selection  rules  are  different  in  general  from  the  linear 
spectroscopy  selection  rules,  TPA  data  may  provide  new  information 
about  electronic  wave  functions  and  energy  levels  in  the  material.  In 
particular,  additional  information  on  band  structure  in  solids  may  be 
obtained.  Practically,  the  TPA  process  has  been  used  to  excite  bulk 
carriers  in  large  semiconductor  crystals  and  generate  color  centers  in 
some  alkali-halides.  These  effects  have  been  used  for  pumping  semi- 
conductor lasers^  and  distributed  feedback  tunable  i,  r.  lasers.  The 
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11.  12 


multi-photon  absorption  provides  a damage  mechanism,  ’ and  thus 
it  may  determine  the  damage  threshold  of  the  material.  This  mechanism 


is  especially  important  in  u.  v.  windows.  With  the  recent  development  in 
high  power  u.  v.  lasers,  practical  interest  has  arisen  in  tlie  absolute  TPA 
coefficients  of  the  window  materials  in  the  u.  v.  region  where  published 


data  are  quite  limited.  In  harmonic  generating  crystals,  the  limited 
conversion  efficiencies  have  been  assigned  to  nonlinear  absorption, 
or  to  phase  mismatch^'^  by  different  authors.  Accurate  measurement  of 
the  nonlinear  absorption  coefficients  is  useful  in  resolving  this  issue. 

In  Sec.  n,  a brief  review  of  pertinent  theories  of  TPA  in  solids  is 
presented.  The  experimental  setup  and  procedure  are  discussed  in 
Sec.  III.  In  Sec.  IV,  the  results  are  presented  and  in  Sec.  V,  our  results 
are  compared  with  calculations  based  on  the  Keldysh  theory  of  high 
frequency  tunneling. 
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ll. Theory 

In  materials  with  band  gap  greater  than  one-photon  but  smaller 
than  the  two-photon  energy  of  an  interacting  light,  it  is  possible  to  have 
absorption  by  simultaneously  absorbing  two  or  more  photons  from  a 
single  laser  beam.  The  two-photon  transition  rate  per  unit  volume  may 
be  obtained  by  the  second  order  perturbation  theory 


2 2^2 
n^c  h 


2ui)  , 


(1) 


where 


2 2 
m u) 


<pfi  ■ 

OD.  - U) 

Ig 
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n is  the  index  of  refraction,  c is  the  velocity  of  light  in  vacuum,  m 

and  e are  the  mass  and  charge  of  the  electron,  respectively,  and  I is 

the  intensity  of  the  laser.  Summations  are  extended  over  final  states 

|f  >,  and  intermediate  states  | i> . The  initial  ground  state  is  denoted  by 

|g>.  a is  a unit  vector  in  the  direction  of  the  electric  field,  p.^  is  the 

matrix  element  of  the  momentum  operator  between  the  states  i and  j. 

hii)--  is  the  energy  difference  between  the  latter  states. 

^ J 

In  the  standard  perturbation  approach,  the  matrix  elements  and 

the  energy  levels  [see  Eq.  (1)]  are  found  from  the  solutions  of  the 

unperturbed  Hamiltonian.  Another  approach  to  attack  the  high  order 

15 

transition  rates  problem  was  developed  by  Keldysh,  who  takes  into 
account  the  main  effect  of  the  electromagnetic  field  at  the  first  step  of 
the  calculations.  In  this  formalism,  the  transition  rate  between  the 


Stark- shifted  states 


i|r(r,t)  = U (r)  exp 
p(t) 
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(2) 


p(t)  • r 


1 «f(p(T))  dT 

•'0 


where 

-«  eS 

p(t)  = P + — sin  cDt  , 

is  obtained.  This  approach  is  useful  since  it  takes  into  account  all  orders 
of  the  perturbation  expansion.  Using  this  technique,  the  multi-photon 

15 

ionization  rate  in  both  gases  and  solids  had  been  calculated  by  Keldysh, 

and  others.  Keldysh's  expression  for  the  transition  rate  reduces  to 

tunneling  at  low  frequencies  while  at  high  frequencies  the  result  reduces 

to  multi -photon  absorption.  In  a solid  with  fundamental  energy  band  gap 

15 

E , using  a simple  parabolic  energy  band  structure,  the  multi-photon 

6 

transition  rate  per  unit  volume  is  given  by  the  Keldysh  formula: 


where 

2 2/  * Z 

A = E + e E /4m  (u  . 
g 

Z i Z + 

$(z)  = e 3 dy  is  the  Dawson  integral.  E is  the  peak  field,  m 

0 

is  the  effective  reduced  mass  of  the  electron  and  hole  pair,  ( ) means 
the  integer  part  of  the  argument. 

The  complex  field  dependence  of  this  equation  may  be  simplified  in 
the  high  frequency  limit.  The  dominant  field  dependence  is  then  given  by 
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the  last  factor  in  the  transition  rate  [Eq.  (3)].  The  N-photon  absorption 
coefficient  is  defined  by  the  expression 


(4) 


which  has  only  a very  weak  field  dependence  through  A.  The  cross 
sections  may  be  obtained  by  transmission  versus  laser  intensity 
measurements.  Xn  the  case  of  one-  and  two-photon  absorbing  medium 
one  obtains,  with  "Yj  = ol  and  ~ transmitted  intensi  y 

1(7, (,t)  = (1  . 

(5) 


where  f is  the  length  of  the  sample,  and  R is  the  reflectivity. 

Multiple  reflections  are  neglected  in  the  derivation  of  the  above  equation. 

Since  we  wish  to  compare,  in  this  paper,  multi-photon  transitions 
of  different  order,  we  have  used  the  Keldysh  theory  for  estimations  of 
multi-photon  cross  sections.  Because  this  theory  is  based  on  a very 
simplified  band  structure,  the  present  calculations  are  not  expected  to  be 
accurate.  Satisfactory  agreement,  however,  was  achieved  between  the 
dispersion  of  P measured  here  and  the  Keldysh's  expression,  using  one 
adjustable  parameter. 

The  TPA  coefficient  may  also  be  expressed  in  terms  of  the 
imaginary  part  of  the  third  order  nonlinear  susceptibility  tensor 
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tu,  u;  , -(u).  For  cubic  crystals,  with  433^  43m  or  m3m  symmetry  and 

( 3 ) 18 

a single  laser  frequency,  x'  only  three  independent  elements. 

These  elements  may  be  obtained  by  the  inea? uremenc  of  the  TPA  cross 

section  in  three  different  polarization  conditions.  First,  it  is  convenient 

to  measure  the  diagonal  element  according  to  the  relation 


P = (32  TT^  m/n^c^)  [3x^j^(-io,  -lu)]  , 


(6) 


which  is  derived  in  the  Appendix.  Here  tu  is  the  laser  angular  frequency, 

c is  the  speed  of  light  in  vacuum,  n^  is  the  index  of  refraction.  We  are 

19 

using  the  Maker  and  Terhune  convention,  with  a degeneracy  factor  3, 
because  there  are  three  distinct  permutations  of  the  positive  and  negative 
frequency  arguments.  As  shown  in  the  Appendix,  for  linearly  polarized 
light  in  the  [1,  1,0]  direction,  the  TPA  coefficient  is 


2 2 2 

p = (IbTT  (u/uqC  )[3x^'j^(-uJ,tu,  u), -uj)  + 6x^'jj.yy(-u),  uj,  cu, -uj) 


+ 3 x"  (-tu,  UJ,  UJ,  -uj)] 

^xyyx  > > > 


(7) 


And  for  circularly  polarized  light,  with  light  propagating  along  [0,  0,  1]  , 


2 2 2 

p = (IGtt  u)/nQC  )[3x]^^(-u),  (D,  oj, -uj)  + 6 x]^'jj.yy(-uj,u),  lu, -tu) 


- 3 x"  (-U),  u),  u),  -co)  ] 

'^xyyx  > > > 'J 


(8) 
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III.  Experimental  Setup  and  Procedure 

The  Nd:YAG  laser  oscillator-amplifier  system  is  described  in 
Ref.  17.  Passive  mode-locking  is  employed  with  a sparx-triggered 
electro- optical  shutter  to  select  a single  pulse  having  energy  up  to  10  mJ . 
The  beam  profile  was  measured  carefully.  It  is  Gaussian  both  in  space 
and  time.  The  pulses  have  an  average  duration  of  30  ps  at  full  width 
half  maximum.  The  fundamental  laser  frequency  is  doubled  by  a 
temperature  tuned  90°  phase  matched  CDA  crystal  (see  Fig.  la).  The 
energy  conversion  efficiency  at  100  MW  input  is  50%.  The  third  harmonic 
(355  nm)  is  generated  by  mixing  the  fundamental  and  the  second  harmonic 
frequencies  in  an  angle  tuned  KDP  crystal.  The  energy  profile  of  the 
beam  is  monitored  by  photo-diode  array  (AR),  and  fitted  with  a Gaussian 
curve  with  a l/e  radius  d = 0.  78  mm.  The  fourth  harmonic  (266  nm) 
is  generated  by  frequency  doubling  the  second  harmonic  in  a temperature 
tuned  90°  phase  matched  ADP  crystal.  The  beam  is  again  fitted  with  a 
Gaussian  of  radius  d = 0.68  mm.  The  overall  energy  conversion 
efficiency  to  355  nm  is  10%,  to  266  nm  15%.  Unwanted  frequencies  are 
filtered  out  by  color  filters  as  shown  in  Fig.  la,  or  by  a dispersive 
prism  used  at  minimum  deviation  in  order  not  to  distort  the  beam. 

In  order  to  determine  the  TPA  cross  section,  the  energy  trans- 
mission coefficient  through  the  sample  is  measured.  For  a Gaussian 
beam  with  a maximum  on-axis  input  intensity  1^,  we  may  write  for  the 
transmitted  intensity 
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I^(l-R)^  exp[-(r/d)  - (t/r)  ]exp(-a^) 

l+pIp(l-R)  exp[-(r/d)^  - (t/T)^][l-exp(-af)] 

while  the  energy  transmission  coefficient  is  given  by: 


(9) 


T = 


out 


in 


CD  « 

Ztt  I r dr  I I(r,  £,  t)  dt 

*'o  -*  


I n' 
P 


d T 


= 2 


ca 

a exp(-a<)  (1-R) | h + £ J (1 -R)[l -exp(-a()) 

[l-exp(-a«l  ■’o  I “ P 


■ exp(-x  ) ( dx  . (10) 


At  low  intensity,  we  find: 


lim  T = (1-R)^  , 

I ^0 
P 


lim  d ” P ^ 

I -*0  P 

P 

The  approximation  is  good  for  1^(1-^  )P/q  < 1 and  a«  1, 

which  is  satisfied  in  all  samples.  Tl  S:  j slope  of  1/T  versus 

at  low  intensities  will  give  the  TPA  coefficient.  During  the  TPA  measure- 
ment a planar  vacuum  photo-diode  Dj  is  used  with  a fast  oscilloscope 
to  make  sure  that  the  beam  is  a single  pulse.  The  fluctuation  in  the  pulse 
area  and  duration  is  taken  into  account  by  monitoring  the  energies  of  each 
pulse  at  1.  06  pm  and  its  harmonics.  A fundamental  Gaussian  pulse 
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I {r,t)  = I exp[-(r/d)^-(t/T)^]  , (13) 

yj 

generates  a second  harmonic 

= gj^Ip  axp[-2(r/d)^-2(t/T)^)  . 


Here  g-.  is  the  conversion  efficiency.  The  energy  of  the  fundamental 
‘='20) 


IS  “I  T A where  A = jr  d . The  second  harmonic  energy  is 

0 p 


1/2  ^ 

= V ' IpT 

7I^/^I^(tA)-,  where  (tA),  = tA  • Similarly,  the  third 
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and  fourth  harmonics  have  ^4u)  " ®4a)  ^ 

- 3/2  - 1 

with  (tA),  = tA  • 3 , (tA)  = tA  • 8 . For  constant  conversion 
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3 1 /2  / 

However,  as  is  shown  in  Fig.  2,  the  experimental  ratio  of  ^ ^ 

depends  somewhat  on  the  input  intensity.  This  deviation  from 
a horizontal  line  indicates  that  some  saturation  effects  in  the  harmonics 
conversion  efficiency  and/or  a change  in  the  laser  pulse  shape  occurs  as 
the  laser  intensity  is  enhanced  by  increasing  the  voltage  on  the  laser 
amplifiers.  Fora  constant  voltage,  however,  the  ratio  (<^o'^‘^3(u^  ^ 
is  nearly  constant  while  the  intensity  is  \*aried  with  attenuators.  The 
latter  effect  and  the  narrowly  distributed  data  shown  in  Fig.  2 show  that 
the  main  fluctuation  is  coming  from  the  fundamental  beam  and  that  TA 
technique  provides  a good  correction  from  pulse  to  pulse  at  the  same 
nominal  intensity.  In  Fig.  3,  we  show  a set  of  data  plotted  with  and 


INTENSITY  (ARBITRARY  UNITS) 


Figure  3. 


Plot  of  the  reciprocal  energy  transmission  coefficient  vs.  intensity 
(dots,  analyzed  with  tA  technique)  or  energy  (cross,  analyzed 
without  tA).  Note  the  two  sets  of  scales. 
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without  a correction  by  the  tA  technique.  This  data  clearly  demonstrate 
the  above  statement. 

The  tA  correction  technique  may  introduce  systematic  error  at 
high  amplifier  voltage  setting  due  to  saturation.  In  the  intensity  range 
of  the  TPA  experiment  performed  here,  we  observed  that  the  average  tA 
increases  by  50%  at  the  peak  power  used  (100  MW)  using  (tA)^  obtained 
from  the  beam  at  2u)  for  intensity  corrections.  We  note  that  the  harmonic 
signal  used  for  obtaining  the  various  tA  were  obtained  from  the  same 
harmonic  beams  utilized  in  the  TPA  experiments.  The  generation  of  3u) 
which  is  done  in  an  angle  tuned  crystal,  critically  depends  on  the 
preceding  second  harmonic  generation  step  which  is  done  in  a 90  CDA 
crystal.  Because  of  the  nature  of  the  mixing  process,  it  is  clear  that  the 
value  of  tA  obtained  from  the  energy  of  the  beam  at  2u)  is  more  reliable 
than  the  value  obtained  from  the  beam  at  3u).  Therefore,  we  analyzed 
our  data  using  (tA)^  for  the  correction. 

The  laser  energy  is  measured  with  the  diode  which  is  calibrated 

with  respect  to  a thermopile  and  a pyroelectric  energy  meter.  The  outputs 
of  all  photo-diodes  were  recorded  by  oscillograms  or  by  homemade  peak 
detect-amplifier-holding  circuits  which  give  the  reading  on  digital 
displays.  For  analyzing  the  TPA  experimental  data,  the  approximate 
value  of  P is  found  using  Eq.  ^12)  and  then  a set  of  curves  are  generated 
using  Eq.  (10)  with  different  values  of  p close  to  the  approximate  value. 

The  TPA  coefficient  p is  then  chosen  by  the  best  fit  curve  as  represented 
by  the  line  in  Fig.  4a. 

Because  the  TPA  cross  sections  are  rather  small,  a lens  (see  Fig.  la) 
is  used  to  increase  the  effective  nonlinear  absorption  by  increasing  the 
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Figure  4.  a.  Plot  of  the  reciprocal  energy  transmission  vs.  maximum  incident 
intensity  for  NaBr.  Sample  thickness  is  0.  74  cm,  laser  electric 
field  is  in  the  [1,0,0]  direction.  The  drawn  theoretical  line  is 
I generated  using  Eq,  (10)  with  p = Z.  5 x 10'^  cm/MW. 

I b.  Typical  plot  of  inverse  energy  transmission  vs.  the  input  intensity 

1 for  materials  with  band  gap  greater  than  ZAui.  Note  the  differences 

; in  the  horizontal  and  vertical  scales  in  parts  a and  b. 
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intensity  of  the  beam  at  th3  sample.  Since  the  geometry  is  critical  in 
the  determination  of  |3,  the  focal  length  of  this  lens  (LI)  is  measured 
carefully.  This  is  done  at  eav.h  wavelength  used  because  the  dispersion 
is  quite  large  in  the  spectral  region  of  the  measurement. 

The  samples  are  put  behind  the  focal  point  of  the  lens  (Ll)  to  avoid 
damage  caused  by  the  effect  of  self-focusing.  Reflections  from  the  back 
surface  of  the  beam  splitter  BS  3 (Fig.  la)  and  the  front  surface  of  BS  5 
are  used  to  avoid  error  introduced  by  TPA  in  the  beam  splitters  them- 
selves. 

In  order  to  eliminate  any  possible  accumulated  effects  on  the  TPA 
measurement  (such  as  F-center  generation),  each  transmission  data 
point  was  performed  at  a fresh  site  in  the  sample.  This  was  done  by 
translating  the  sample  after  each  laser  shot. 

The  circularly  polarized  beam  used  to  probe  the  anisotropy  of 

io  generated  by  passing  the  beam  through  a Fresnel  rhomb,  with 
proper  beam  polarization.  The  circularity  is  checked  by  a Gian  prism 
analyser. 

The  uncertainty  of  the  relative  ^alue  of  p measured  here  is  quite 
small.  The  main  contribution  to  this  uncertainty  is  the  uncertainty  in 
the  experimental  value  of  tA.  On  the  average  this  introduces  an 
uncertainty  of  ± 15%.  The  scattering  of  the  data  points  gives  on  the 
average  another  ± 10%.  The  individual  uncertainty  for  each  sample  is 
lited  in  Tables  1 and  2.  The  systematic  error  in  the  absolute  determination 
of  p , an  uncertainty  which  should  be  added  to  the  previous  one,  is  ± 7% 
in  absolute  energy  calibration,  ± 28%  in  the  average  value  of  tA  at  low 
intensity,  and  ± 8%  caused  by  the  ±1.5  mm  uncertainty  in  the 
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determination  of  the  focal  length  of  the  lens  (Ll).  Thus,  the  absolute 
uncertainty  sums  up  to  ± 43%. 

The  experimental  setup  for  the  relative  measurement  with  the  two- 
channel  normalization  technique^  is  shown  in  Fig.  lb.  A double  Fresnel 
rhomb  is  used  here  to  vary  in  an  efficient  way  the  intensities  going  to  the 
sample  and  reference  channels.  The  ratio  of  the  energy  transmitted 
through  the  sample  and  the  reference  is  given  by 
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At  proper  filter  factor,  F = integrands  in  Eq.  (15)  are  equal 

and  the  ratio  p will  be  independent  of  intensity.  V/ e have  then 
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(16) 


To  check  the  reliability  of  the  setup  for  this  experiment  we  used  identical 
samples  of  KJ.  as  sample  and  reference.  Fig.  5 shows  clearly  that  the 
expected  value  of  F = 1 was  obtained. 
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IV.  Experimental  Results 

At  355  run,  we  measured  the  absolute  TPA  coefficients  in  three 
alkali-halides,  as  well  as  in  KDP  and  10  other  harmonic  generating 
crystals.  These  are  all  listed  in  Table  1.  Upper  limits  for  p in  the 
commonly  used  u.  v.  window  materials.  Sapphire  and  fused  silica,  with 
band  gaps  greater  than  2Au)  (7,  02  eV)  are  also  tabulated.  For  the 
materials  measured  with  band  gaps  greater  than  2fix  , the  transmission 
T is  practically  constant  to  intensity  up  to  20  GW/cm  . Fig.  4b  shows 
a typical  transmission  plot  for  material  with  large  band  gap.  Higher 
beam  intensities  could  not  be  used  due  to  the  onset  of  surface  damage. 

The  anisotropy  of  the  TPA  was  measured  in  Rbl.  Linear  polarized 
light  in  the  [1,,0,0]  and  [l,f,0]  directions  show  the  same  value  of  p. 
According  to  Eqs.  (6)  and  (7),  this  suggests  the  isotropy  condition 

2y"  + y"  = X"  . (17) 

■^xxyy  ^yyx  '^xxxx 


The  circularly  polarized  beam  showed  a smaller  TPA,  The  ratio  for  the 

(3) 

three  elements  of  x'  calculated  from  Eqs.  (6-8)  are  xlL™  • x'l-  • 

- 1 3 3 

x"  = 1:0.455:0.09,  with  the  absolute  value  x"  = 2.  72  X 10  cm  /erg. 
^xyyx  ’ ^xxxx  ' ^ 

Anisotropy  of  the  tetragonal  crystal  CD  A was  also  observed.  The  value 
of  p with  the  beam  polarized  in  the  x-y  plane  is  larger,  by  a factor  1.  8, 
than  that  obtained  when  the  E-field  is  parallel  to  z axis.  Although  some 
second  harmonic  generation  may  occur  in  the  former  case,  this  radiation 
is  not  phase  matched.  The  power  conversion  efficiency  at  an  input 
intensity  of  10  GW/crn  is  estimated  to  be  10  and  therefore  the 
measured  value  of  p is  not  affc  ed. 
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table  2 


Material 

Eg(eV)' 

Surface^ 
Normal  N 

E/I 

e| 

c 

P^^(KW)“ 

t(cm) 

0(  cm/MW  )'^ 

* / ® 
m /m 

a 

n 

NaCi 

6.36 

(1  0 0) 

10 

0 

11 

16 

0.74 

3.5  X 10“^  t 25t 

5.3  X lo"^ 

1.637 

NaCi 

6.36 

11  0 0) 

10 

1 

11 

16 

0.74 

3.6  X lo"^  i 30\ 

NaBr 

7.7 

(1  0 OJ 

10 

0 

11 

8 

0.74 

2.5  X lo"^  t 15t 

1.03  X 10”^ 

1.803 

KCi 

6.52 

11  0 01 

10 

0 

11 

20 

0.74 

1 .7  X lo”^  ± 20t 

2.46  X 10"^ 

1.571 

KCi 

11  0 01 

10 

0 

11 

20 

0.74 

1 .8  X lo"^  + 30t 

(F-centers) 

KCi 

6.52 

11  0 01 

10 

0 

11 

20 

2.10 

2.7  * lO"^  * 30t 

KBr 

6.00 

11  0 01 

10 

0 

11 

4.5 

0.74 

2.0  X lo”^  ± 30t 

1.46  X 10”^ 

1.76 

KI 

5.06 

11  0 01 

10 

0 

11 

6 

0.47 

3.75  X 10‘^  ± 30t 

4.8  X 10'^ 

1.938 

RbCt 

7.29 

|1  0 Jl 

10 

0 

11 

0.6P 

1.02  X lo"^  i 15t 

4.54  X 10”^ 

1.576 

RbCt 

7.29 

11  0 01 

10 

0 

11 

1.91 

1.26  X 10“^  i 20\ 

RbBr 

5.40 

11  0 01 

10 

0 

11 

1.0 

2.18  X 10“^  ± 20t 

2.02  X lo"^ 

1.686 

Rbl 

5.04 

11  0 01 

10 

0 

11 

1.0 

2.49  X lo’^  ± 20t 

1.15  X lo”^ 

1.917 

CaCOj 

5.ee 

0.43 

2.4  X 10“’  ± 30t 

1.11 

1.750  dip) 

1.526 

(n  ! 

e 

KDP 

6.95 

1 .061J->0. 53U 

67 

0.57 

2.7  X 10"“*  t 30% 

0.72 

1.560  (iip) 

1.510 

<"e' 

cut 

ADP 

6.81 

0.5 

2.4  X lo"'’  * 30% 

0.90 

1.580  (rip) 

1.527 

(n 

e 

AljOa 

7.3 

C-axis 

0.5 

2.7  X lo"*  ± 30% 

-4 

0.353 

1.833 

(no) 

At,0, 

7.3 

1 C-axis 

0.5 

2.7  X 10  ±30% 

2 3 

-5 

SiO, 

7.8 

70 

0.64 

4.5  X 10 

50 

(fused  silica) 

SiOj 

7.8 

“ C-axis 

1.2 

4.5  X lo"^ 

50 

(puarti) 

LiF 

11.6 

(1  0 01 

(0 

0 

11 

73 

0.47 

<2  X 10 

c 

1.415 

CaFj 

10.0 

(1  0 01 

(0 

0 

11 

60 

0.47 

<2  X 10 

Table  1.  The  TPA  experimental  results  at  355  nm. 

a.  Energy  gand  gap  (E  ),  self-focusing  critical  power  {Pcr)> 
and  refractive  inde:f  (n)  are  given  whenever  available. 

b.  All  measurements  were  done  normal  to  the  crystal  surface. 

c.  £/l£|  is  the  direction  of  the  laser  polarisation. 

d.  The  total  uncertainty  in  p is  the  sum  of  the  relative 
uncertainty  (listed  in  the  table  for  each  material)  and  the 
systematic  uncertainty  of  the  experiment  estimated  to  be 
± 43%. 
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In  addition  to  the  absolute  determination  of  (3  we  measured  the 
relative  values  of  (3  in  RbBr  versus  KI  and  Rbl  versus  KI  as 
described  above.  The  ratios  obtained  p(RbBr )/p(KI)  = 0.  305  and 
p(RbI)/(3(KI)  = 0.  668  agree  within  10%  with  the  absolute  values  in  Table  1. 
This  agreement  between  the  two  measurements  suggests  that  the 
uncertainty  estimation  presented  in  Sec.  Ill  is  quite  conservative. 

At  266  nm,  we  measured  the  TPA  in  8 alkali-halides  with  band  gap 
less  than  2^a'  (9.  36  eV),  and  in  KDP,  ADP,  Calcite,  Sapphire  and  Si02. 
The  window  materials  LiF  and  CaF^  in  which  no  TPA  is  expected  are 
also  measured.  The  measured  values  are  listed  in  Table  2.  For  alkali- 
halides  the  p values  measured  are  of  the  order  of  10  cm/MW,  while 
for  KDP,  ADP,  calcite  and  sapphire,  they  are  of  the  order  of  10  cm/MW 
For  the  materials  with  band  gaps  greater  than  Zh^jD,  an  upper  limit  for  p 
is  given.  It  is  important  to  note  that  fused  silica  and  quartz  have  quite 
small  TPA  at  this  wavelength. 

Fluorescence  is  observed  in  all  the  alkali-halides  along  the  path  of 
the  beam.  Coloring  can  be  seen  with  the  bare  eye  after  a few  shots. 

In  KCl  the  absorption  of  a monitoring  beam  at  532  nm  by  the  color  track 
is  observed  and  can  be  explai  ' \ by  the  spectral  position  of  the  absorption 
band  of  the  F-center.  The  physical  process  involves  the  generation  of 
electron-hole  pairs  by  two-photon  excitation.  The  pairs  recombine  to 
form  excitons  with  the  emission  of  fluorescent  light.  The  exciton  e.iergy 
is  given  up  to  form  and  F-centers.  We  have  found  that  the 

F-centers,  generated  by  several  hundred  shots  of  the  laser  at  the  same 
site  of  the  crystal,  did  not  influence  the  transmission  of  the  crystal  at 
the  laser  frequency  after  they  had  relaxed  in  few  seconds. 


-28- 


table  1 


Material 

E^(eV)  * 

Surface_^, 
Nonnal  N 

E/IeI 

£(cni) 

e (cm/MW) 

d 

a 

n 

KI 

5.06 

11 

0 

01 

(0  0 1)  11 

0.47 

7.29 

X 

10-' 

± 

20\ 

1.765 

RbBr 

5.40 

[1 

0 

01 

(0  0 1) 

1.0 

2.43 

X 

0 

1 

* 

20% 

1.600 

Rbl 

5.04 

(1 

0 

0] 

(0  0 11 

1.0 

5.08 

X 

»-• 

0 

1 

+ 

15% 

1.731 

Rbl 

5.04 

(1 

0 

01 

(0  1 1) 

1.0 

5 08 

X 

lo"' 

4 

2 5% 

Rbl 

5.04 

(1 

0 

01 

Circularly 

polarized 

1.0 

4.62 

X 

1 

0 

1 

20% 

ADA 

45“ 

2 

cut 

2-axiB 

T 6 

3.53 

X 

10-' 

4 

30% 

ADP 

6.81 

45» 

2 

cut 

z-axiE 

2.0 

6.8 

X 

10-' 

4 

35% 

1.548 

<"o> 

1.499 

(n  ) 
e 

CDA 

45* 

2 

cut 

Xz-axiB 

1.4 

8.02 

X 

10-' 

4 

30% 

1.60 

'"o' 

1.57 

'"e’ 

CD*  A 

45' 

2 

cut 

z-aucis 

1.4 

5.12 

X 

10-' 

4 

30% 

1.59 

'"o' 

1 57 

'"e’ 

CD*  A 

45' 

2 

cut 

z-axiB 

1.15 

2.81 

X 

10-' 

4 

35% 

KDA 

45' 

2 

cut 

z-axiB 

1.0 

4.84 

X 

10*' 

4 

30% 

1.59 

("o' 

1.54 

'"e’ 

KD*A 

45' 

2 

cut 

z-axis 

1.0 

2.66 

X 

10-' 

f 

35% 

KDP 

6.95 

45' 

2 

cut 

Z-axiB 

1.0 

5.9 

X 

-0-' 

4 

35% 

1.531 

'"o' 

1.486 

'"e’ 

KD*P 

45' 

2 

cut 

Z-axiB 

1.0 

5.4 

X 

10-' 

i 

35% 

1.53 

("o' 

1.49 

'"e' 

RDA 

45' 

2 

cut 

z-axis 

1.25 

4.99 

X 

10-' 

4 

30% 

1.60 

'"o' 

1.55 

'"e' 

RDP 

45' 

2 

cut 

z-axis 

1.0 

<5.9 

X 

10-' 

1.53 

'"o’ 

1.50 

'"e' 

A 2O3 

7.3 

C-axi: 

5 

0.5 

<1.6 

X 

10-' 

1 796 

("o’ 

SiO 

7.8 

0.64 

<1.25 

X 

10-' 

1.476 

(fused  sil 

ica) 

Table  2.  The  TPA  experimental  results  at  266  nm. 

a.  Energy  band  gap  (Eg),  self-focusing  critical  power  (P^^), 
and  refractive  index  (n)  are  given  whenever  available. 

b.  All  measuvements  were  done  normal  to  the  crystal  surface. 

c.  £/|£l  is  the  direction  of  the  laser  polarization. 

d.  The  total  or  ertainty  in  p is  the  sum  of  the  relative 
uncertainty  (listed  in  the  table  for  each  material)  and  the 
systematic  uncertainty  of  the  experiment  estimated  to  be 
± 43%. 

e.  m*/m  is  the  ratio  of  the  effective  mass  fitted  with  the 
Keldysh  formula  to  the  observed  value  of  P and  the  mass 
of  the  free  electron. 
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Since  the  power  used  in  the  TPA  experiment  is  much  higher  tha,n 

1 Z 26 

the  critical  power  for  self-focusing,  ’ it  is  important  to  look  for  any 

change  in  beam  size  due  to  self-focusing.  The  checks  have  been  done 

with  different  sample  thickness,  also  by  observing  the  stable  color  track 

formed  in  a U-center  doped  KCi  under  microscope.  An  attempt  was 

made  to  record  a single  shot  beam  profile  in  three  dimensions  by  the 

27 

TPA  generated  color  centers  in  KCI,  but  the  coloring  effect  is  not 

strong  enough  to  be  observ'ed  after  a single  shot.  No  evidence  of  self- 

focusing  has  been  observ'ed  so  far,  in  spite  of  predictions  based  on 

26 

critical  pow'er  v'alues  for  self-focusing  extrapolated  by  the  scaling 
_ 2 

factor  0)'  from  X = 1.06  jam  measurements.  Since  the  measurements 

were  done  vdth  Ztm  close  to  the  energy  gap,  the  latter  extrapolation 

might  be  incorrect  because  of  the  positive  or  negative  two-photon  resonant 

contributions  to  the  nonlinear  refractive  index  n^.  In  addition,  the  TPA 

effect  which  consumes  energy  from  the  peak  intensity  location  of  the  beam 

competes  \vith  the  self-focusing  process. 

Our  values  of  P are  about  a factor  5 smaller  than  those  obtained 

28  29 

more  than  a decade  ago  by  Frohlich  and  Park.  ’ With  the  less  well 
controlled  beams,  nominal  values  for  the  intensity  are  probably 
considerably  lower  than  the  actual  intensity  in  hot  spots  due  to  spatial 
and  temporal  fluctuations.  This  could  explain  their  systematically  high 
value  of  p.  A more  recent  measurement  of  TPA  coefficient 
p = 2.  36  X 10  cm/MW,  by  Mollenauer  et  al.  , in  KCI  is  more  in  line 
v/ith  our  data.  Reintjes"^^  has  very  recently  reported  very  low  values  of 
the  TPA  coefficients  in  ADP  and  deuterated  KD  P.  His  value  for  KD  P 
P = (2.  7 ± 0.  7)  X 10"^  cm/MW  is  an  order  of  magnitude  lower  than  our 

_4 

value  for  the  undeuterated  KDP.  His  value  for  ADP  P = (1.  1 ± 0.3)  x 10 
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cm/MW  is  a factor  of  2.  2 lower  than  ours.  The  large  deviation  of  the 
results  in  the  two  latter  crystals  may  be  explained  by  differences  in 
impurity  concentration  in  the  sample  used.  It  is  interest'ng  to  note  that 
for  the  alkali-halides,  we  measured  at  both  u.  v.  frequencies,  the  TPA 
coefficients  at  355  nm  are  of  the  same  order  but  somewhat  larger  than 
those  at  266  nm. 
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V . Comparison  with  the  Keldysh  Theory 

The  magnitudes  for  4-photon,  3-photon,  and  2-photon  absorption 
coefficients  of  NaCl  derived  from  the  Keldysh  theory  [Eq.  (4)]  are 
plotted  as  a function  of  photon  energy  in  Fig.  6.  For  a fixed  N, 
goes  through  a very  shallow  maximum  as  a function  of  photon  energy. 

The  maximum  occurs  at  an  energy  just  before  the  next  lower  order 
photon  process  becomes  possible.  The  band  gap  dependence  and  the 
effective  mass  dependence  of  which  represent  the  material  dependence 

in  the  Keldysh  t>-  -lory,  are  shown  in  Fig.  7. 

The  observed  values  of  (3  in  Table  2 with  the  known  band  gap  for 
each  material  can  all  be  fitted  with  the  Keldysh  theory,  by  adjusting  m 
in  each  case.  These  adjusted  values  of  m are  listed  in  Table  2.  It 
should  not  be  expected  that  these  values  correspond  to  those  determined 
by  transport  properties  from  carriers  near  the  bottom  of  the  conduction 
band.  In  fact,  m in  the  Keldysh  theory  is  largely  a measure  for  the 
average  strength  of  the  dipole  matrix  element. 

Using  the  m*  found  at  266  nm  for  KI,  RbBr  and  Rbl,  the  frequency 
dependence  of  P in  each  of  these  materials  was  calculated.  The  results 
are  plotted  in  Fig.  8,  along  with  the  experimental  points.  It  is  interesting 
to  see  that  the  Keldysh  theory  does  predict  a decrease  of  p with 
increasing  photon  energy  in  qualitative  agreement  with  the  experimental 
trend. 

The  effect  of  multi-photon  absorption  on  breakdown  can  also  be 

calculated  by  the  Keldysh  formula.  In  order  to  compare  this  model  with 

12  31 

the  experimental  data  observed  by  Smith  et  al.  , ’ we  calculated  the 

threshold  electric  field  needed  to  generate  enough  free  electrons  by 


(MKS  UNIT) 


PHOTON  ENERGY  f\(v  (eV) 


6.  Log  of  the  multi- photon  absorption  coefficients  in  NaCl  as 
given  by  the  Keldysh  formula  vs.  photon  energy  (m*/m  = 1 
was  used).  The  steps  indicate  transitions  from  4-  to  3-  and  to 
2-photon  process.  The  arrows  indicate  the  photon  energy  as 
fraction  of  the  energy  band  gap,  Eg. 
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multi-photon  absorption  alone,  and  melt  the  sample  by  plasma  heating. 
Assuming  no  recombinations  during  a square  picosecond  laser  pulse,  the 
temperature  rise  by  the  free  electron  absorption  is  given  by 


at 


2 

e T 
4mc 


W E (At)‘ 


1.22 

1+U)  T 


(18) 


Here  c is  the  specific  heat,  m is  the  free  electron  mass,  T is  the 
momentum  transfer  collision  time,  At  is  the  pulse  duration,  E is  the 
peak  electric  field  strength,  u)  is  the  angular  frequency  of  the  light,  and 
W is  given  by  the  Keldysh  formula  [Eq.  (3)].  Using  Eq.  (18)  with  m 
obtained  from  the  theoretical  fit  of  the  TPA  measurements  at  266  nm 
and  a relaxation  time  t = 5 x 10  sec,  the  electric  field  threshold 
for  breakdown  is  calculated.  The  results  for  KDF  and  fused  si-ica  are 
presented  along  with  the  available  experimental  data  in  Fig.  9. 

As  discussed  in  Ret.  12,  the  breakdown  threshold  at  low  frequency 
is  determined  by  an  avalanche  multiplication  process.  In  this  case  a 
monotonic  increase  in  the  breakdown  threshold  is  expected  with 
increasing  laser  frequency.  This  is  reflected  by  the  experimental  data 
on  fused  silica,  but  the  contribution  to  the  breakdown  from  the  3 -photon 
absorption  process  at  the  highest  frequency  is  appreciable. 

In  KDP,  the  measured  breakdown  thresholds  at  1.  06  pm  and  532  nm 

are  almost  equal.  An  increase  by  a factor  of  two  occurs  at  355  nm.  The 

Keldysh  theory  gives  breakdown  threshold  values  close  to  the  experimental 

data  at  1.06  pm  and  532  nm.  In  this  region  the  relative  contributions  of 

the  avalanche  ionization  and  the  multi-photon  process  to  the  carriers 

••  I-  I-  1 32 

generation  have  to  be  considered  as  discussed  by  Braunlich  et  al. 


THRESHOLD  ELECTRIC  FIELD  (MV/cm) 


PHOTON  ENERGY  (eV) 


Figure  9.  Theoretical  calculation  according  to  Eq.  (17)  (open  points)  and 
experimental  data  (filled  points)  of  the  threshold  peak  electric 
field  for  breakdown  in  KDP  (squares)  and  fused  silica  (triangles). 
The  dashc  curves  are  used  to  aid  in  reading  the  figure.  The 
breakdown  Hi  silica  is  apparently  determined  by  avalanche 
ionization,  and  not  by  multiphoton  absorption. 
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The  larger  experimental  threshold  observed  at  355  nm  might  reflect  the 
nonlinear  loss  in  the  sample,  which  was  not  taken  into  account  in  the 
determination  of  the  breakdown  threshold  reported  in  Ref.  31.  For 
example,  with  TPA  coefficient  p = 6 x 10  ^ cm/MW  and  tightly  focused 
input  beam  with  intensity  I = 2 x 10  W/cm  , the  percentage  loss  per 
100  pm  is  as  high  as  12%.  Therefore,  nonlinear  losses  which  are  not 
negligible  at  the  high  laser  intensity  used  should  be  carefully  considered. 

In  fused  silica,  the  expected  nonlinear  loss  at  355  nm  (3-photon 

4 

absorption)  is  quite  smell.  Therefore  the  threshold  measured  is  reliable. 

When  the  damage  threshold  is  determined  mainly  by  the  avalanche 
mechanism,  multi-photon  absorption  before  the  focal  region,  where  break- 
down would  occur  is  an  advantage.  It  lov'ers  the  intensity  at  hot  spots, 
therefore,  the  possibility  of  self-focusing  and  damage  is  reduced. 

In  conclusion,  we  have  measured  the  absolute  TPA  coefficients  of 
several  commonly  used  u.  v.  materials  and  harmonic  generating  crystals 
at  two  frequencies.  The  frequency  dependence  is  adequately  described  by 
the  theory  of  Keldysh.  TPA  has  no  appreciable  effect  on  the  conversion 

efficiency  in  KDP  and  its  isom^rphs.  We  concur  in  this  conclusion  reached 

14 

by  Reintjes  et  al. 
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APPENDIX 


The  relations  between  the  TPA  coefficient  in  different  polarization 

(3) 

conditions,  and  the  corresponding  components  of  x (-lu,  cu,  ou, -cu)  are 
obtained  by  defining  the  third-order  nonlinear  polarization.  In  cubic 

(3) 

materials  only  4 components  of  y independent.  By  using  the 

(3) 

general  permutation  rules  for  y the  additional  symmetry  because 

we  are  dealing  with  only  a single  laser  frequency,  we  get  for  the  folloviung 
three  cases: 


A.  the  electric  field  parallel  to  the  crystal  x axis  with  amplitude  E^, 


the  electric  fi  'd  in  the  [1,  1,0)  direction  with  same  amplitude  E^, 

f3  E'""  + 6 x^^^  E E e"'  + 3x^^^  e’" 

■^xxxx  X X '^xxyy  y x y '^xyyx  y x 

3 x^^^  E*  + 6 xl^^  E E e’%  3x^^^  E^  E* 

'^xxxx  y y ■'xxyy  x y x ^xyyx  x y 


C.  the  electric  field  is  circularly  polarized,  = Eq(x+  i y)  , 
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(A3) 


I p • £ ^ I 

The  TPA  work  done  on  the  material  is  given  by  w = Re  j — ^ 
TPA  coefficient  is  defined  by 


P = wVl^  . 


(A4) 


Using  the  expression  for  the  laser  intensity  for  linearly  polarized  light 


I - ^ 

^ ■ 877  0 ’ 

and  for  circularly  polarized  light 


(A5) 


2 


(A6) 


with  the  nonlinear  polarizations  given  above,  Eqs.  (6-8)  given  in  the  text 
are  easily  obtained. 
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C,  SUMMARY  OF  WORK  ON  SUPERCONDUCTING  MATERIALS. 

Professor  M.  Tinkham,  Principal  Investigator. 

Three  entirely  different  classes  of  novel  superconducting  materials 

1-3 

have  been  studied,  and  publications  have  appeared. 

The  first  of  these  are  the  superconducting  layered  compounds,  such 
as  TaS-  and  NdSe-,  both  with  and  without  organic  intercalates  such  as 
pyridine.  These  compounds  have  highly  anisotropic  superconducting 
properties,  the  critical  field  parallel  to  the  plane  of  the  layers  being 
orders  of  magnitude  higher  than  the  perpendicular  critical  field.  There 
appeared  to  be  some  potential  for  utilizing  the  high  parallel  critical  field 
with  possibly  strong  pinning  of  flux  lines  between  the  layers  to  make 
superconductors  with  unique  high-field  properties.  One  of  the  critical 
questions  to  be  resolved  was  that  of  dimensionality;  does  the  material 
act  like  a collection  of  independent  2-dimensional  superconducting  planes, 
or  is  the  interplane  coupling  strong  enough  to  yield  3 -dimensional,  albeit 
highly  anisotropic,  behavior.  Prober’s  experiments  were  aimed  at 
answering  such  questions  by  measuring  the  critical  field  as  a function  of 
angle  with  the  plane  and  as  a function  of  temperature,  and  also  of 
measuring  the  diamagnetism  above  T^  associated  with  thermodynamic 
fluctuations.  The  latter  appeared  to  offer  particularly  direct  information 
about  dimensionality,  since  the  temperature  dependence  in  the  2-D  and 
3-D  cases  are  different.  In  fact,  although  the  measurements  were  more 
extensive  and  careful  than  any  previous  ones,  sample-related  problems 
set  limits  to  the  definiteness  of  the  conclusions.  In  particular,  warping 
and  crinkling  of  the  planes  made  the  notion  of  precise  parallelism 
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unrealistic.  Nonetheless,  from  the  data  we  could  conclude  that  the 
materials  do  act  3-dimensional  over  a considerable  temperature  range 
on  either  side  of  T^;  in  fact,  no  convincing  evidence  for  2-dimensional 
behavior  was  found.  Also,  the  very  weak  dependence  of  interplane 
coupUng  on  the  thickness  of  the  organic  intercalating  layer  shows  that 
the  coupling  is  not  of  a simple  tunneling  type,  which  would  have  an 
exponential  variation.  A paper  on  this  work  by  D.  E.  Prober  (now  at 

Yale)  and  M.  R.  Beasley  (now  at  Stanford)  has  been  published. 

2 

The  work  of  Powell,  Skocpol  and  Tinkham,  was  designed  to  pursue 
the  feasibility  of  making  superconducting  materials  of  superior  pinning 
strength  by  tailoring  the  scale  of  the  microstructure  to  the  intervortex 
spacing,  typically  hundreds  of  angstroms.  The  materials  studied  were 
NbN  and  NbC  and  intermediate  compositions.  Fine  particles  of  a Nb- 
containing  compound  were  prepared  by  spray  drying  a suitable  solution. 
These  particles  were  then  reduced  to  the  desired  compound  by  heating 
them  in  an  atmosphere  of  niixed  with  CH^  or  NH^.  They  were  sub- 
sequently pressed  and  sintered  into  solid  rods,  whose  superconducting 
properties  (T  , H _,  J ) were  measured.  A wide  range  of  conditions 
was  explored,  varying  the  spray  drying  technique,  solution  concentration, 
etc.  , varying  the  temperature  and  gas  mixture  for  the  heat  treatment, 
and  then  varying  the  final  sintering  conditions.  With  a certain  amount  of 
optimisation  of  these  many  variables,  it  proved  possible  to  make  NbN 
with  pianing  properties  (as  measured  by  J^)  somewhat  superior  to  the 
best  reported  on  samples  made  by  more  conventional  techniques. 
However,  it  was  hard  to  retain  optimum  grain  size  while  obtaining 
sufficient  sintering  of  particles  to  make  a cohrent  solid.  Thus,  although 
this  work  basically  supports  our  initial  theoretical  notions,  this  technique 
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does  not  appear  to  offer  much  promise  for  the  fabrication  of  practical 
superconductive  materials.  The  fine  particle  preparation  technique 
may  have  wider  applicability,  however. 

The  final  class  of  materials  worked  on  is  that  of  superconductive 
mic ro- composite  wire,  produced  metallurgically  in  the  way  first  reported 
by  Tsuei.  This  wire  is  typically  made  by  melting  together  Nb,  Sn,  and 
Cu,  quenching  at  such  a rate  that  small  Nb  grains  precipitate  from  the 
Cu-Sn  solution.  Upon  drawing,  the  Nb  grains  become  extremely  long, 
thin  filaments  (of  finite  length)  embedded  in  a Cu-Sn  solid  solution  matrix. 

A suitable  heat  treatment  allows  the  Sn  to  react  with  the  Nb  filaments, 
forming  a sheath  of  the  high  performance  superconductor  Nb^Sn.  This 
yields  a wire  containing  thousands  of  filaments  of  micron  diameter  in  the 
cross  section  of  the  wire.  Because  of  the  fineness  of  the  filaments,  the 
composite  is  ductile,  despite  the  hardness  of  the  Nb^Sn.  Also,  being 
made  by  purely  macroscopic  metallurgical  methods,  it  should  be  relatively 
inexpensive,  and  of  interest  for  applications,  such  as  fusion  confinement 
magnets,  requiring  large  amounts  of  superconductive  wire.  An  obvious 
question,  however,  was  whether  such  a wire  containing  finite  length 
superconductive  filaments  a copper  matrix  was  truly  superconductive. 

On  the  other  hand,  a simple  tl.euretical  model  showed  that  even  if  the 
copper  were  fully  normal,  the  resistance  would  be  reduced  by  a factor  of 
order  10^  below  that  of  the  copper  matrix,  just  because  of  the  presence  of 
the  very  long -filamentary  superconductive  segments.  The  expected 
residual  resistance  would  be  below  the  level  detectable  by  conventional 
instrumentation,  and  in  fact  low  enough  that  it  would  not  interfere  with 
practical  application.  Davidson  developed  a superconductive  voltmeter 
sensitive  enough  to  detect  the  expected  residual  resistance.  What  he  then 
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found  was  that  at  low  currents,  the  material  is  a true  superconductor, 
the  proximity-effect  superconductivity  in  the  copper  being  sufficient  to 
carry  the  current  from  one  filament  to  the  next.  At  higher  currents 
(and/or  magnetic  fields)  he  found  a very  small,  but  finite  resistance, 
similar  to  that  expected  theoretically.  Moreover,  he  pointed  out  the 
need  to  invoke  notions  from  percolation  theory  to  obtain  an  appropriate 
understanding  of  the  observations.  That  is,  above  a percolation 
threshold  estimated  to  be  15%  concentration  of  superconductor  by  volume, 
one  expects  a continuous  superconducting  path  by  mere  random  close 
contacts  between  filaments.  This  work  is  at  present  continued  with 
support  from  other  sources. 
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D.  SUMMARY  OF  WORK  ON  FRACTURE  MECHANICS. 

Professor  J.  W,  Hutchinson,  Principal  Investigator. 

Work  was  completed  in  two  areas  during  the  funding  period: 

(i)  large  scale  yielding  crack  problems  and  (ii)  near  tip  elastic-plastic 
fields  for  growing  cracks. 

Extensive  numerical  results  were  generated  for  fully  plastic  plane 
stress  crack  problems  for  two  of  the  most  important  crack  configurations, 
the  center  cracked  strip  and  the  edge  cracked  strip  in  bending.  These 
solutions  can  be  used  in  an  approximate  way  to  construct  simple  formulas 
governing  fracture  conditions  for  the  full  elastic -plastic  configuration. 

The  formulas  reduce  at  low  load  to  the  elastic  results  and  at  high  load 
they  merge  into  the  fully  plastic  solution.  To  assess  the  accuracy  of  the 
simple  formulas  some  full  scale  finite  element  calculations  were  carried 
out.  The  agreement  with  the  simple  formula  predictions  was  excellent. 

2 

This  work  is  reported  in  the  paper  by  C.  F.  Shih  and  J.  W.  Hutchinson. 

For  the  first  time,  stress  and  strain  fields  have  been  found  at  the 
tip  of  a crack  growing  steadily  in  an  elastic-plastic  strain  hardening 
material.  The  material  is  characterized  by  a flow  theory  of 
plasticity  together  with  a bilinear  effective  stress -strain  curve.  All 
three  cases  of  interest  were  considered,  i.  e.  , anti-plane  shear,  plane 
stress  and  plane  strain.  Perhaps  the  most  important  finding  is  that  a 
small  amount  of  strain  hardening  leads  to  a relatively  strong  singularity 
compared  to  the  ideally  plastic  case.  It  is  still  true  that  the  strain  fields 
at  the  crack  tip  are  not  nearly  as  singular  in  the  growing  case  as  in  the 
stationary  case.  But  the  theoretical  results  suggest  that  it  still  may  be 
possible  to  employ  the  concept  of  a growing  crack-tip  stress  intensity 
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factor  to  analyze  stably  growing  cracks  in  analog  to  the  stress  intensity 

concept  that  has  been  so  successful  in  treating  stationary  cracks.  This 

work  is  summarized  in  the  paper  by  J . C.  Amazigc  and  J.  W. 

3 

Hutchinson. 
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